Soliton-radiation trapping in gas-filled photonic crystal fibers 
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We propose an optical trapping technique in which a fundamental soliton traps an ultrashort 
small amplitude radiation in a symmetric hollow-core photonic crystal fiber filled with a noble gas, 
preventing its dispersion. The system is Raman- and plasma- free. Trapping is due to the cross 
phase modulation effect between the two pulses. The trapped radiation inside the soliton-induced 
potential will oscillate periodically due to the shock effect, similar to the motion of a mechanical 
pendulum. 

PACS numbers: 42.65.Tg, 42.81.Dp 



Introduction - Optical trapping has been observed 
in a variety of different situations since the invention of 
low-loss optical fibers [lj. Self-trapping is at the core 
of the existence of optical solitons, where the nonlinear 
self-phase modulation (SPM) is used to compensate the 
group velocity dispersion (GVD) Bright and dark 
solitons have been observed in the anomalous and the 
normal dispersion regimes, respectively, in a number of 
elegant experiments [3|-|5(- Trapping can also take place 
between a pair of pulses, when they lie in opposite dis- 
persion regimes to satisfy the group velocity matching 
condition. For instance, red-shifted solitons due to intra- 
pulse Raman can trap across the zero dispersion wave- 
length weak blue-shifted dispersive wave-packets on the 
femtosecond scale In an unconventional situation, 

a bright soliton can also be trapped in the normal disper- 
sion, when it is coupled to a dark soliton in the anomalous 
dispersion regime via the nonlinear cross-phase modula- 
tion (XPM) effect [H,[l4|. 

The invention of photonic crystal fibers (PCFs) is a 
milestone in the field of nonlinear fiber optics [ljl [lij . 
Hollow-core (HC) PCFs have been exploited to explore 
the nonlinear interaction of light in a new wide-range 
of distinct media. HC-PCFs based on a Kagome lat- 
tice have been successfully used recently in the investi- 
gation of light-matter interactions in presence of gases 
171 ]. Filling these fibers by a Raman-active g dras- 
tically reduction in Raman threshold has been obtained 
fl8| . High harmonic, and efficient deep UV generation 
have been also demonstrated in HC-PCFs filled by noble 
(Raman- inactive) gases [ljl, [2(J. Using intense broad- 
band pulses to access the ionization regime of the no- 
ble gases, soliton self-frequency blueshift has been ob- 
served and thoroughly studied in these fibers 2lM24l] . On 
the other hand, launching intense narrowband pulses can 
lead to ionization-induced asymmetric SPM and univer- 
sal modulational instability |25j . 

In this Letter, we show that optical trapping between 
two different short pulses can exist, when both pulses 
propagate together in the anomalous dispersion regime 
of a HC-PCF filled by a noble gas. The intensities of the 



two pulses are below the ionization limit of the gas. 

Principle of operation — Consider the propagation of 
two short pulses in a guided nonlinear-Kerr medium. As- 
suming that only one of the two pulses is intense enough 
to induce nonlinear phenomena, and the dispersion is 
anomalous for both pulses. Hence, the strong pulse can 
maintain its temporal shape due to the induced SPM, 
while the other pulse will suffer from strong dispersion 
due to the absence of compensation between nonlinear 
and dispersive effects. Only if the two pulses can prop- 
agate with the same group velocity, the nonlinearity of 
the strong pulse will affect the second pulse via XPM 
that can compensate the dispersion-induced broadening, 
and the weak pulse will propagate as a temporally lo- 
calized radiation. In analogy with quantum mechanics, 
the strong pulse creates a potential well that traps the 
weak pulse inside it, and both pulses travel together with- 
out suffering from dispersion-induced broadening [lll.ll2[. 
If the nonlinear medium is Raman- or plasma-free, the 
dispersion must be anomalous for both pulses to allow 
dispersion-compensation to take place. 

System design - To observe this temporal trap- 
ping due to the soliton-induced potential, we require a 
nonlinear-Kerr medium, where two different ultrashort 
pulses lying in the same anomalous dispersion regime can 
acquire a common group velocity. In previous works, the 
group velocity matching condition was achieved when a 
red-shifting soliton was in the anomalous dispersion and 
the small amplitude wave was in the normal dispersion 



regime [llj, [12| . One possible technique for our proposal 



is to use a solid-core PCF with multiple zero dispersion 
wavelengths. However, this type of fibers requires high 
fabrication tolerance. Moreover, undesirable emission of 
resonant radiation will be out of control in these fibers. 
Another technique is to launch the two pulses in two dif- 
ferent polarization states in an elliptical solid core fibers. 
However, relatively long pulses must be used in this case 
to minimize the effect of Raman - and it is that case one 
lacks the attracting feature of trapping ultrashort pulses. 

A HC-PCF filled by a noble (Raman- inactive) gas and 
with a symmetric core is considered to be the perfect host 
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to demonstrate our proposal. This can be achieved by 
launching two different pulses with the same frequency, 
however in two different circular polarization states, in 
the deep anomalous dispersion regime of the fiber. Hav- 
ing the same frequency allows the two pulses to have the 
same group velocity, as well as to be naturally in the same 
dispersion regime. Using different polarization states al- 
lows to monitor the two pulses at the fiber output. The 
reason for using circular polarization states is to avoid 
unwanted additional nonlinear phenomena to emerge if 
linearly polarized states are injected 0]. 

Governing equations for soliton-radiation trapping - 
To study the induced trapping between a propagating 
strong pulse and a small amplitude radiation with the 
same central frequency and different circular polarization 
states in a lossless nonlinear medium, the following set of 
coupled nonlinear Schrodinger equations can be used Q , 
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where z is the longitudinal coordinate along the fiber, 
t is the time coordinate, Aj (z,t) is the slowly varying 
amplitude of the jth pulse, Dj(id t ) = ifiijd t — ^/3 2 jd 2 , 
fiij and /?2j are the first- and second-order dispersion 
coefficients of the jth pulse (j = 1,2), higher order dis- 
persion coefficients are neglected, 7 is the nonlinear Kerr 
coefficient, and A/3 is the fiber birefringence. A\ and 
A 2 correspond to the strong pulse and the weak radia- 
tion, respectively. In our proposed system, the medium 
is Raman- and plasma-free, f3±i = j3\2, P21 = /?22, and 
A/3 = for perfectly symmetric cores. In each equa- 
tion, the first (second) nonlinear term represents the self- 
(cross-) phase modulation effect. We introduce the fol- 
lowing rescalings: £ = z/zo, t = t/to, ipi = A\/Aq, 
4> 2 = A 2 /A , Al = 3/(2 7Z0), z a = t\j |/3 2 i|, where t is 
the input pulse duration. In a reference frame moving 
with a group velocity v g , which is defined as v g = l//3n, 
the two coupled equations can be replaced by 
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where the nonlincarity of the weak radiation is ne- 
glected. The solution is well-known for the first equa- 
tion as a nonlinear Schrodinger soliton, ib\ (£, r) = 
77 sech (jit) exp (ir7 2 £/2) , with amplitude r) [2|. We as- 
sume that the solution of the second equation has the 
form ^(^t) = ipo f (f) exp (iq£), where ipo and q are 
the amplitude and the propagation constant of a pulse 
with temporal profile / (t) . Substituting this form in 
the second equation, we obtain a linear one-dimensional 
Schrodinger problem in time, 



d 2 T f + U(T)f = -qf, 



(3) 



tial well, / (r) is the cigenfunction, and —q is the cor- 
responding eigenvalue that represents a discrete energy 
level. We have found the analytical expressions of the 
modes of this potential well. The fundamental (even) 
mode is /o (r) = sech Po (ar) with 5=5 V 2 Poi and po = 
\ [— l + ^/Tf\. The first-order (odd) mode behaves as 
/i(t) = sech Pl (ar)tanh(aT), with q = \ rj 2 (2 - 3pi), 
and p\ — po — 1. For this system, there are only two local- 
ized modes. In principle, ipo for either the fundamental or 
the first-order mode can take an arbitrary value as long 
as ipo -C 77, since these modes propagate linearly. The 
problem can also be solved numerically by using sparse 
matrices (26[. The obtained numerical results validate 
the analytical formulas for the eigenfunctions and eigen- 
values given above. It is worth mentioning that an at- 
tempted solution of a similar problem is presented in an 
eminent quantum mechanics textbook [271 ] . however, a 
full analytical form of the eigenfunctions is presented for 
the first time in this work. 

Propagation of a dispersive radiation in an argon-filled 
HC-PCF is depicted in Fig. Q] Panels (a,b) show the 
dispersion-induced broadening of the radiation in absence 
of the soliton- induced potential, when the launched pulse 
takes the shape of either the fundamental or the first- 
order mode, respectively. In presence of a fundamental 
soliton [with a temporal and spectral evolution shown in 
panels (c,d)], an XPM- induced trapping for the radiation 
can take place. Panels (e,f) and (g,h) present the prop- 
agation of a localized radiation, when it takes the form 
of the fundamental or the first-order mode, respectively. 
Higher-order dispersion coefficients and self-steeping ef- 
fects are neglected in this simulation. 

Shock effect — Since the nonlinear medium is Raman 
and plasma free, and the two pulses are supposed to be 
arbitrarily short, the self-steeping (shock effect) will have 
a significant role in the pulse dynamics. In this case, Eq. 
([2|) becomes 

id^ipi + \d 2 ipi + \ipi\ 2 ipi = -iT sh d T ^i\ 2 ipij , 

i&^2 + \ 8^2 + 2 |Vl I 2 V>2 = ~i 2 T sh 8 T (|^l| 2 & 

'(4) 

where r s h = (w to) 1 is the normalized shock coeffi- 
cient, and wo is the input pulse central frequency. In 
order to study the shock effect on both the soliton and 
the radiation, we will adapt the variational perturba- 
tion method for our problem Q. Using this method, 
the RHS of Eqs. (01 are considered as perturbations for 
the previously obtained solutions of ipi and ip2. Hence, 
the soliton and the radiation will maintain their func- 
tional shape during propagation, however, their param- 
eters will vary. For both pulses, one has to obtain the 
Lagrangian C = /_ Cd dr, where Cd is the Lagrangian 
density. Using the reduced Eulcr-Lagrange equations, a 
set of ordinary differential equations that tracks the spa- 
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FIG. 1. (Color online). (a,b) Propagation of a dispersive ra- 
diation solely in an argon-filled HC-PCF, when the radiation 
takes the form of (a) the fundamental mode fo (t), or (b) the 
first-order-mode /i (r). (c-h) Propagation of a fundamental 
soliton and a dispersive radiation together in an argon-filled 
HC-PCF. The group-velocity matching condition is satisfied. 
The two pulses have the same central wavelength A = 800 
nm, and opposite circular polarizations. (c,d) Temporal and 
spectral evolution of the soliton. (e-h) Temporal and spectral 
evolution of the radiation, when it is in the fundamental mode 
(e,f), or in the first-order-mode (g,h). Other simulation pa- 
rameters are r\ = 8.9, tpo = 0.5, and to = 50 fs. Higher-order 
dispersion coefficients and self-steepening are neglected. 



tial evolution of the pulse parameters can be obtained. 
Since this technique is well illustrated in the literature 
for the soliton case Q, we will show only its final re- 
sult, Vi(£i T ) = ??sech [rj (r — /3'£)] exp [irj 2 ^/2), with 
P' = T s h i] 2 ■ This means that the soliton acquires only 
a linear temporal shift during propagation due to the 
shock term. For the radiation, 



-lm{0 2 d e ^+2e^} 



2|^i| 2 |^ 2 | 2 , 
(5) 

where e = — 2r s h d T (jV'il 2 ^2) , Im stands for the imagi- 
nary part. Assuming that the radiation is in the funda- 
mental mode, 

i>2 (£, r) = -00 scch Po [a (r - r p )] exp [-id (r - r p )] , (6) 

where t p and 5 are the shifts due to the applied pertur- 
bation in the temporal and spectral position of the pulse 
center, respectively After calculating the Lagrangian 
and applying the reduced Eulcr-Lagrange equations, one 



can show that the following equations describe the vari- 
ation of the parameters: 
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where I = 0Fr(p o ) /V (p + |), .9 (r) = 
Po tanh [a (r — t p )], T is the gamma function, and 
Re stands for the real part. These coupled equations 
can be solved numerically. 

The effect of the shock term on the radiation is de- 
picted in Fig. [2] which shows the temporal and the spec- 
tral evolution of the dispersive radiation. As shown in 
panel (a), the radiation acquires a varying group veloc- 
ity, which differs from that of the soliton. The reason 
is related to the factor 2 that distinguishes between the 
SPM and XPM processes [see Eq. Q]. This group veloc- 
ity difference between the two pulses causes the observed 
periodic oscillations of the radiation inside the soliton. 
The dynamics of the radiation is similar to the motion 
of a mechanical pendulum, or an oscillating particle in- 
side a potential. Initially, the particle is settled at the 
potential minimum. Due to the shock effect, it acquires 
a velocity that is highest in correspondence with that 
minimum. As the 'particle' (which is in our case the 
small amplitude radiation) moves upward, it undergoes 
a deceleration, resulting in a spectral blueshift as shown 
in panel (b). When the particle reaches its maximum 
position, it stops and changes its direction downwards, 
acquiring an acceleration that corresponds to a spectral 
rcdshift. Passing by the potential minimum, the parti- 
cle goes upward to the other side of the potential. This 
periodic motion continues during the propagation along 
the fiber. In the case of the fundamental mode, there 
is a good qualitative agreement between the analytical 
and numerical results, see Fig. [2{a,b). The odd mode 
behaves similarly to the even mode as portrayed in pan- 
els (c,d). However, the analytical results fail to mimic 
the numerical ones in this case. We attribute this to the 
reason that the radiation does not preserve its temporal 
shape during propagation. The two-lobes that constitute 
the radiation packet exchange energy, similar to partial 
coupling between adjacent waveguides, and this prevents 
the perturbation theory to be a good approximation to 
the dynamics. 

Conclusions - We present a trapping technique be- 
tween a strong fundamental soliton and a weak disper- 
sive radiation in a symmetric HC-PCF filled by a noble 
gas. The medium is Raman- and plasma-free. Trapping 
is due to the cross-phase modulation effect between the 
two pulses having the same central frequency but oppo- 
site circular polarization states. Although both pulses lie 
in the deep anomalous dispersion regime of the fiber, this 
configuration satisfies unconventionally the group veloc- 
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FIG. 2. (Color online). Propagation of a fundamental soliton 
and a dispersive radiation in an argon-filled HC-PCF in pres- 
ence of self-steepening effect. Other simulation parameters 
are similar to those used in Fig. [T] Temporal and spectral 
evolution of the radiation, when it is in the fundamental mode 
(a,b), or in the first-order mode (c,d), respectively. Solid black 
lines are the numerical solutions of the coupled differential 
equations Eq. ([7|). 



ity matching condition. From quantum mechanics point- 
of-view, the soliton creates a potential well that traps the 
particle-like radiation inside, and the two pulses travel 
together without suffering dispersion-induced broaden- 
ing. This potential well allows only for two bound states. 
Exact analytical formulas for the fundamental and the 
first-order modes of the radiation have been obtained. 
The soliton shock effect leads to a periodic oscillation 
of the radiation inside the potential, similar to a pen- 
dulum motion. As a consequence, the radiation central 
frequency is shifted periodically. 
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